function Figure = fun_createFigure(F, A, B, B1, freq, I_row, I_col, ...
            noise_intensity, noise_mean_num, gamma,...
    phase_shift, phase_shift_tilts, filepath)

        
    
    % 生成均匀相移的矩阵
    Figure = cell(F,1);
    Delta = cell(F,1);
    for f = 1 : F
        Delta{f,1} = phase_shift(f)*ones(I_row,I_col);
    end
    % 生成非均匀相移的矩阵
    Alpha = cell(F,1);
    Delta_mac = Delta;Q = 1;
    for f = 1 : F
        Alpha{f,1} = [0, f, f, f, f, f] * phase_shift_tilts;
        [~,Delta_mac{f,1}] = Maclaurin(I_row,I_col,Q,Alpha{f,1});
    end


    % 均匀相移+非均匀相移 = 最终加在物体上的相移
    for f = 1 : F
        Delta{f,1} = Delta{f,1} + Delta_mac{f,1};
        Delta{1,1}(:) = 0;
        addDelta = Delta{f,1}(670:760,270:315);
        save([filepath,'-',int2str(f),'-addDelta.mat'],'addDelta');
        mean_phase_shift = sum(addDelta(:))/(I_row*I_col);
        save([filepath,'-',int2str(f),'-mean_phase_shift.mat'],'mean_phase_shift');
    end

  
    for f = 1:F
        delta = Delta{f,1};
        for c = 1:I_col
              Phi = c*2*pi*freq/I_col.*ones(I_row,1);
              phi_delta = Phi+mod((f-1),4)*pi/2+delta(:,c);
              Figure{f,1}(:,c) = A + B*cos(phi_delta)+B1*cos(2*phi_delta);
        end
    end

%     num = 1;
    for f = 1 : F
%        str = [filepath,'-',int2str(f),'.bmp'];
       Figure{f,1} = Figure{f,1}.^(1/gamma);
       Figure{f,1} = imnoise(Figure{f,1},'gaussian',noise_mean_num,noise_intensity);
%        imwrite(Figure{f,1},str);fprintf('------%s条纹图生成！------ \n', str);
%        num = num + 1;
    end
end